Formation of the Pd atomic chain in hydrogen atmosphere 
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Abstract 

The formation of a Pd atomic chain in a hydrogen atmosphere was investigated by measurements 
of conductance and vibrational spectroscopy of a single molecular junction, and the theoretical 
calculation. While atomic chains were not formed for clean 3d and 4d metals, in the case of Pd (a 
4d metal) atomic chains could be formed in the presence of hydrogen. Stable atomic chains with 
two different atomic configurations were formed when the Pd atomic contact was stretched in a 
H2 atmosphere; highly conductive short hydrogen adsorbed atomic chain and low conductive long 
hydrogen incorporated atomic chain. 

PACS numbers: PACS numbers: 73.63.Rt, 73.40.Cg, 73.40. Jn 
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INTRODUCTION 



Chains of single metal atoms forming bridges between metal electrodes have attracted 
wide attention due to their unique properties |l( and their potential applications as ultimate 
conducting wires in nanoscale electronic devices 2|. In the case of Au, atomic chains up to 
2 nm long have been created in UHV, and their physical properties have been investigated 
with various techniques (2 - 4 ] . The formation of atomic chains was also reported for other 
5d metals, such as Pt and Ir. On the other hand, it has not been possible to prepare atomic 

u n 

chains in the case of 3d and 4d metals for the following reasons |2J, |5[. The formation of 
the atomic chain requires relatively strong metal-metal bonds within the chain compared 
to other bonds in regions close to the chain. In the case of 5d metals, the bonds in the 
atomic chain, in which all atoms are at the surface, is strengthened due to the relativistic 
effects of the valence electrons . In contrast, the stabilization of the bond strength in 
atomic chain made from 3d and 4d metals is not significant, making the fabrication of chains 
from these metals extremely difficult. Despite these difficulties, however, various interesting 
properties have been predicted. For example, theoretical studies have shown that the Pd 
atomic chain is ferromagnetic, even if the length of the atomic chain is short (~3 atom) 
Consequently, successful fabrication of atomic chains of 3d and 4d metals, Pd in particular, 
is of great scientific interest. 

Recently, it was shown that atomic or molecular adsorption on metal surfaces can stabilize 
the metal atomic contact. 2 nm long Ag atomic chains were created in the presence of oxygen 
in UHV at 4 K, while clean Ag only forms short chains 16] . In the case of Pd, Pd atomic 

n 

chains were stabilized under the hydrogen evolution reaction in solution [7j. In UHV, the 
conductance of the Pd atomic contact decreased from 2Gq (Go = 2e 2 /h to 0.5Go due to 
the introduction of hydrogen jq]. While it was revealed that it is possible for the metal 
atomic chain to be stabilized in the presence of other atoms and molecules, the structure 
of the stabilized atomic chain is not clear. It is not known whether a molecule adsorbs on 
the chain surface or if it is incorporated into the chain. Vibrational spectroscopy of a single 
molecular junction is one promising technique used to investigate the atomic configuration 
of the single molecular junction [9|, |lOj. In the present study, we report on the formation 



of Pd atomic chains in a hydrogen atmosphere. The structure of the atomic chain has 
been verified by comparing conductance measurements and vibrational spectra of a single 
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molecular junction to those from theoretical calculations. 



EXPERIMENTAL 

All measurements were performed using a mechanically controllable break junction 
(MCBJ) technique (see Ref.[2| for a detailed description). A notched Pd (0.1 mm diam- 
eter, 99.99% purity) wire was glued onto a bending beam and mounted in a three-point 
bending configuration inside a vacuum chamber. Once under vacuum and cooled to 4.2 K, 
the Pd wire was broken by mechanical bending of the substrate. The bending was then 
relaxed to form atomic-sized contacts between the wire ends using a piezo element for fine 
adjustment. H2 was injected into the chamber via a capillary. DC two-point voltage-biased 
conductance measurements were performed by applying a voltage in the range from 10 to 
150 mV. AC voltage bias conductance measurements were performed using a standard lock- 
in technique. The conductance was recorded for a fixed contact configuration using a 1 mV 
AC modulation at a frequency of 7.777 kHz while slowly ramping the DC bias between -100 
and +100 mV. 

EXPERIMENAL RESULTS 

Figure [1] (a) shows typical conductance traces for clean Pd contacts and for Pd contacts 
in a H2 atmosphere, recorded with a bias voltage of 0.1 V. For the clean Pd contacts, the 
conductance decreased in a stepwise fashion and the corresponding conductance histogram 
(Fig. QJb)) showed a peak near 2 Go, which corresponds to a clean Pd atomic contact. After 
the introduction of H 2 , a plateau below 1 Go was frequently observed and the corresponding 
histogram (Fig. [1] (b)) showed a peak around 0.4 Gq. The obtained conductance histogram 
agreed with the previously reported results BJ. Here it should be noticed that the Pd atomic 
contact in a H2 atmosphere could be stretched to quite long lengths (~0.5 nm), which 
suggested the formation of an atomic chain. In order to investigate the chain formation, 
the length histogram was investigated for the last plateau. Figure [2]^a) shows the length 
histogram of the last conductance plateau for the Pd contacts in a H 2 atmosphere, together 
with that for the clean Pd contacts. The length of the last plateau was taken as the distance 
between the points at which the conductance dropped below 1.0 Go and 0.05 Go for the Pd 
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contact in a H2 atmosphere, and the boundaries were 2.5 Go and 1.0 Go for the clean Pd 
contact. The length histogram for the Pd contacts in a H2 atmosphere did not change when 
the upper boundary was changed from 1.0 Go to 2.5 Go- The length histogram showed that 
the Pd atomic contact could be stretched as long as 0.5 nm in a H 2 atmosphere, while the 
clean Pd contact broke within 0.3 nm. Since the Pd-Pd distance is 0.27 nm for bulk Pd, the 
final 0.5 nm long plateau corresponds to an atomic chain of approximately two Pd atoms. 

The formation of the Pd atomic chain was supported by the analysis of the conductance 
trace. Fig. [2] (b) shows the average return lengths as a function of chain length, recorded 
with a bias voltage of 0.1 V. The return length was defined as the distance over which the 
two electrodes were moved back to make contact after the junction broke. The return length 
was averaged over 2000 breaking cycles. Apart from an offset of 0.1 nm due to the elastic 
response of the banks, the average return length was approximately proportional to the 
chain length, suggesting that a fragile structure was formed with a length corresponding to 
that of the last plateau {2, [5]. The long plateau over 0.5 nm length and the one by one 
relationship between the chain and return length indicated that a Pd atomic chain could 
be formed in a H2 atmosphere. Figure 09(c) shows the average conductance of the atomic 
chain as a function of the chain length. The curve was obtained by adding all measured 
conductance traces from the start value (1 Go) onward, and dividing each length by the 
number of traces included at that point. It was observed that the mean conductance of the 
atomic chain decreased with increasing chain length. 

In order to investigate the atomic configuration of the atomic chain, the vibrational 
spectroscopy of a single molecular junction was measured for the Pd atomic contacts in a H2 
atmosphere. Figure [3](a) shows examples of the d 2 I/dV 2 curves (vibrational spectroscopy of 
a single molecular junction) for the Pd atomic contacts in a H2 atmosphere taken at a zero 
bias conductance of 0.15~0.6 Go- The spectral shape and energy of the vibrational modes 
varied with the contacts. Peaks and dips were observed in the spectra depending on the 
conductance of the contacts. When the conductance of the contact was lower than 0.4 Go, 
only peaks were observed in the spectra. Conversely, both peaks and dips were observed in 
the spectra when the conductance of the contact was larger than 0.4 Go- 

The peaks in the spectra indicated a conductance enhancement above a certain voltage, 
while the dips indicated a conductance suppression. The conductance enhancement can be 
explained by the opening of an additional tunneling channel for electrons that lose energy to 
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vibrational modes 



121 ] . Likewise, the conductance suppression could be explained as the 



backscattering of electrons that lose energy to vibrational modes in the ballistic contacts in 
which the electron transmission probability is close to one. When a molecule symmetrically 
couples to both metal electrodes and the number of the conduction channel is one, the 
theoretical investigation predicts that the conductance of the contact would be enhanced by 
phonon excitation for contacts having a conductance below 0.5 Go, while the conductance 



would be suppressed for the contacts having a conductance above 0.5 Go [Hi, [12|. The 
experimental result that only peaks were observed for the contacts having conductance 
below 0.4 Go agrees with this theoretical prediction. Furthermore, the results that both 
peaks and dips were observed for contacts with conductances above 0.4 Go, indicate the 
presence of more than one conductance channel. In the latter case, the presence of more 



than one channel makes it difficu 
by the phonon excitation 
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t to discuss the direction of the conductance change induced 



Our observations indicate that the Pd atoms bridge between the Pd electrodes for contacts 
with conductances above 0.4 Go, as discussed below. If a hydrogen bridges the gap between 
the Pd electrodes, electrons would be transmitted through Is atomic orbital (in the case 
of a hydrogen atom bridge) or bonding or anti-bonding molecular orbitals (in the case of 
a hydrogen molecule bridge) which would be modulated by Pd. Since these orbitals are 
not degenerate and are energetically well separated, only one orbital can contribute to 
the electron transport through the hydrogen atom or molecule bridge. Since the present 
experimental results suggested multiple conduction channels, hydrogen can not have formed 
a bridge between the Pd electrodes as electrons would have been transmitted through Pd d 
orbitals. 

In the spectra of Fig. [3X the vibrational modes were observed at 25 meV (conductance of 
contact: 0.15 G ), 35 and 75 meV (0.2 G ), 45 meV (0.4 G ) and 45 meV (0.6 G ). While 
the energy of the vibrational modes varied with the contacts, the distribution shows some 
defining characteristics. Specifically, it depended on the conductance value of the contacts. 
The distribution of the vibrational energy of the contacts having conductance below 0.3 Go 
was different from that above 0.3 Go- Figure 3(b) and (c) show the distribution of vibrational 
energy for contacts having conductances below 0.3 G and 0.3~0.6 G , respectively. The 
histograms were obtained from 250 spectra of different contacts. The histogram showed 
two peaks at 29 meV (14 meV in width) and 64 meV (25 meV in width) for contacts 
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with conductance below 0.3 Go and broad single peak at 35 meV (32 meV in width) for 
contacts with conductance of 0.3~0.6 Go- Since the conductance and vibrational energy 
of the atomic contacts is strongly correlated with their atomic configurations, the present 
results indicate that atomic chains with two different atomic configurations would be formed 
during stretching of the Pd contacts in a H 2 atmosphere. The conductance of the atomic 
chain decreased with the chain length, as shown in Fig. [T](a) and Fig. [2] (c). The poorly 
conductive long atomic chain shows two vibrational modes around 29 and 64 meV, while the 
highly conductive short atomic chain shows one vibrational mode around 35 meV. Here it 
should be noted that the spectra for the contact of 0.2 Go (Fig. Efa)) shows two vibrational 
modes, indicating the presences of two vibrational modes in the contact. The double peaks 
in the histogram did not, however, originate from the two atomic contacts with different 
vibrational modes and similar conductance values. Briefly, we comment on the features of the 
conductance histogram. While the vibrational spectra suggest the formation of atomic chains 
with two different atomic configurations, the conductance histogram shows a single feature 
around 0.4 Go- We think that single feature in the conductance histogram is composed of 
two broad peaks. The formation of atomic chains with two different atomic configurations 
is supported by the conductance trace. The last curve in Fig. UK a) shows steps around 0.5 
Go and 0.2 G . 



CALCULATION RESULTS 



To gain insight regarding the structural and vibrational modes of our Pd chain in a H 2 at- 
mosphere, we performed density functional theory (DFT) calculations using the Gaussian03 
program Q. The hybrid DFT approach, B3LYP Q, was employed with the Stuttgart- 
Koln effective core potentials for Pd [15[ and Dunning's cc-pVDZ basis set for H [I6j. The 
atomic contact part was modeled using two types of clusters, (a) Pd4-Pd-Pd4 and (b) Pd4- 
Pd4, where a hydrogen molecule was introduced in the area around the central Pd atoms. 
The geometry was fully optimized for the respective models and normal mode analyses were 
performed for the obtained structures. Through preliminary calculations, we have found 
that H2 cannot attach to the central part of Pd chains as a molecule, and easily dissociates 
to hydrogen atoms on the Pd cluster. Thus, geometry optimizations were performed for the 
respective model clusters with the adsorption of two hydrogen atoms and we succeeded in 
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locating minimum energy structures with no imaginary frequencies for both model clusters. 

Figure H] shows the vibrational modes containing the movement of hydrogen atoms calcu- 
lated for the optimized structures of (a) Pd4-H-Pd-H-Pd4 (referred to as "short Pd chain") 
and (b) Pd3-Pd-H-H-Pd-Pd3 (referred to as "long Pd chain"). The H 2 molecule dissociates 
to H atoms in both structures {thh ~0.192 nm in structure (b)) which work to stabilize 
the contact part of the Pd chain. The present calculations reveal that the short hydrogen- 
adsorbed shows the 5 vibrational modes of the Pd chain around 23-31 meV, while the long 
hydrogen incorporated Pd contact shows the vibrational modes at 41 and 59 meV. 

The vibrational modes observed in the present experiments are discussed based on these 
calculation results and previously reported results of nanocontacts. Pd nanocontact and 
Pd nanocontact in a H 2 atmosphere were investigated using point contact spectroscopy. In 
these experiments, Pd crystal vibrational modes were observed around 15-20 meV for the 
Pd nanocontact and, likewise, the vibrational modes of dissolved H atoms was observed 
around 60 meV for contacts in a H2 atmosphere. In the present study, the histogram of 
the vibrational modes showed two peaks around 29 meV and 64 meV for the long atomic 
chain, and a broad single peak around 35 meV for the short atomic chain. Based on these 
experimental and theoretical calculation results, an experimentally obtained broad single 
peak around 35 meV for the short atomic chain would correspond to the 5 calculated vibra- 
tional modes ranging from 23-31 meV. Since energy differences among the five modes are 
small, they appear as a single broad peak in the histogram. 35 meV is close to the energy of 
the vibrational mode of the Pd crystal (15-20meV) and the vibrational modes of Pd in the 
stem parts of the electrodes might contribute to the 35 meV peak. The two experimentally 
obtained peaks around 29 meV and 64 meV for the long atomic chains correspond to the 
calculated vibrational modes at 41 meV and 59 meV. For the long atomic chains, the vibra- 
tional modes of Pd in the stem part of the electrodes and the modes of dissolves H atoms 
in the Pd crystal can contribute to the peaks around 29 meV and 64 meV, respectively. 
In both short and long atomic chains, the peaks in the histogram of the vibrational modes 
were broad. The atomic chains whose atomic configuration was slightly modified from the 
stable structure could be formed under the present experimental conditions, and thus, the 
vibrational modes from all of these structures can contribute to the peaks, causing them to 
broaden. Here, we briefly comment on the short atomic chains. As discussed before, the 
vibrational spectra suggest that electrons can be transmitted through the Pd rather than 



7 



the hydrogen for contacts with conductance above 0.4 Go (short Pd chain), which agreed 
with the theoretical calculation results. 

The conductance of the Pd atomic contacts with similar atomic configurations has been 
investigated using first-principle m olecular dynatnies —ens Q. The conductance 
values of short and long Pd chains are reported as 1 Go and 0.5 Go, respectively. It was 
found that the conductance of the short Pd chain is smaller than that of the long Pd chain. 
Our experimental results show atomic chain conductance decreased with chain length, as 
shown in Fig. Ufa) and Fig. EJ^c), which is in agreement with the results of theoretical 
calculations. The smaller conductance values obtained in the experimental results can be 
explained as arising from the adsorption or incorporation of hydrogen onto the stems of the 
Pd electrodes since hydrogen molecules may dissociate and migrate into the bulk Pd. The 
theoretical calculations show that the conductances of the hydrogen-adsorbed and hydrogen- 
incorporated Pd contacts decreased by 30~60% depending on the amount of the hydrogen 
in the Pd electrodes [18]. Our experimental and calculated results combined with previous 
theoretical simulations 17| revealed that the atomic chain of two Pd atoms in contact with 
the hydrogen bridge can be formed when a Pd contact is stretched in a H2 atmosphere while 
clean contacts do not, in fact, form an atomic chain. 



CONCLUSIONS 



We have investigated the atomic configuration and conductance of Pd atomic contacts 
in a H 2 atmosphere. The formation of Pd atomic chains was demonstrates using the length 
histogram of the last conductance plateau. Vibrational spectroscopy and the results of 
theoretical calculations show the formation of highly conductive short hydrogen-adsorbed 
Pd chain and low conductive long hydrogen incorporated Pd contacts. 
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FIG. 1: (a) Breaking and return conductance traces for clean Pd contacts (dotted line), and for Pd 
contacts in a H2 atmosphere (line), (b) Conductance histograms for the clean Pd contacts (dotted 
line), and for Pd contacts in a H2 atmosphere (line). Each conductance histogram was constructed 
from 1000 conductance traces recorded with a bias voltage of 0.1 V during the breaking of the 
contact. 
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FIG. 2: (a) Length histogram for clean Pd contacts (dotted line) and Pd contacts in a H2 atmo- 
sphere (line).(b) Average return lengths as a function of chain length, (c) Average conductance as 
a function of chain length. Each curve in (a)-(c) was constructed from 2000 conductance traces 
recorded with a bias voltage of 0.1 V during the breaking of the contact. 
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FIG. 3: (a) Typical d 2 I/dV 2 curves for Pd contacts in a H2 atmosphere taken at conductance 
ranging from 0.15~0.6 Go. (b) The distribution of vibrational energy for the contacts with con- 
ductance below 0.3 Gq. (c) The distribution of vibrational energy for contacts with conductance 
ranging from 0.3~0.6 Gq. 
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FIG. 4: Vibrational modes containing the movement of hydrogen atoms calculated for the optimized 
structures of (a) Pd 4 -H-Pd-H-Pd 4 (short Pd chain) and (b) Pd 3 -Pd-H-H-Pd-Pd 3 (long Pd chain). 
The length of the arrows indicates the amplitude of the modes. 
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